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Abstract 

In this paper we offer an alternate route for investigating soliton solutions 
in hydrogen-bonded chains. This is done by examining a class of systems of 
two coupled real scalar fields. We show that this route allows investigating 
several models for hydrogen-bonded chains in a unified manner. We also 
show how to investigate interesting issues, in particular the one concerning 
classical or linear stability of solitonic solutions. 
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1 Introduction 



Hydrogen-bonded materials that appear in condensed matter, and in organic 
and biological systems, may transport energy and charge, and this transport 
seems to be due to proton transfer along the hydrogen-bonded chains |], 0, 
13- B 01 • T ne conduction via proton migration along the chain in hydrogen- 
bonded materials appears as follow. We can think of a hydrogen-bonded 
material as being genericaly represented by 

•••X-H---X-H---X-H--- 

where (— ) and (• • •) stand for the covalent and the hydrogen bond, respec- 
tively. Here X represents a group of atoms - for ice it is OH, where O stands 
for oxygen. The hydrogen bond bridges the system in many repetitions of the 
unit cell X — H- • •, and the system is usually considered as a unidimensional 
macroscopic chain. 

In such a macroscopic chain, protons are transfered via two distinct mech- 
anisms. The intrabond mechanism accounts for proton migration through 
the bridge, and this proton migration changes the role of the covalent (— ) 
and the hydrogen bond (•■•)• 111 this case there is a local disturbance of the 
neutral charge of the chain, and so this defect is called an ionic defect. The 
interbond mechanism accounts for the interchange of position of the covalent 
bond. Here there is a local orientational change of the covalent bond, and so 
this defect is called an orientational (or Bjerrum) defect. 

Hydrogen-bonded materials have been studied by many authors. A par- 
ticularly interesting investigation of soliton solutions for hydrogen-bonded 
chains was introduced in Ref. ||. In this case, the longitudinal dynamics 
of protons was investigated within the context of one-component models for 
hydrogen-bonded chains. 

In a more recent work a specific model of hydrogen-bonded chain was 
examined @. Here the model is of the two-component type, and is based 
on the assumption that the coupling between the proton sublattice (the first 
component) and the acoustic mode of the heavy-ion sublattice (the second 
component) is linear. In this case the solitonic solution is a kink or a bell- 
shape soliton, depending of the specific parameters which control the proton 
on-site potential, which is represented by a symmetric (kink) or asymmetric 
(bell-shape) double-well potential. 
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In the present paper, we get inspiration from recent works || [HJ to offer 
an alternate route for examining soliton solutions in hydrogen-bonded chains. 
This will be done by introducing a class of systems of two coupled scalar 
fields. Here, however, we shall work with the same assumption considered 
in 0, which states that the coupling between the proton sublattice and the 
acoustic mode of the heavy-ion sublattice is linear. 

As we are going to show, the class of systems we shall introduce can be 
used to map the model investigated in and other models, in a unified 
manner. Within this context, our investigations provide a new route for 
examining soliton solutions in hydrogen-bonded materials. Before doing this, 
however, let us first introduce the class of systems of two coupled fields. 



2 Two coupled real scalar fields 

The investiagtions we are now going to do is inspired on former works || 
|9|, [HJ, from where we borrow notation, which is standard notation O] in 



Field Theory. The class of systems we shall investigate is described by the 
following Lagrange density 

C = 1(1 + a^d^d^ + h 2 d aX d a X + abd Q <j)d a x - U ((/>). (1) 

We are working in (1 + 1) dimensions. Hence a = 0, 1, and x a = (x° = 
t,x l = x), and x a = (x — t, X\ = —x). Here and % are rea l scalar 
fields, dimensionless, and a and b are real and dimensionless parameters. 
The potential U(<j>) is chosen as 

and is defined from the function H = H(<f)), which is a smooth but otherwise 
arbitrary function of the field (p. The reason for the above choice (0) for the 
potential will appear later. 

In the above model, evidently, each function H(<p) defines a specific sys- 
tem, and so we have a class of systems of two coupled real scalar fields. The 
coupling between the two fields is a derivative coupling - the third term in 
(HD - and is the same for every specific system one introduces by some specific 
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function H(4>). We then recognize that systems belonging to the above class 
are different when they differ on the self-interaction term, or better, on the 
potential for the scalar field <j). 

The Euler-Lagrange equations that follow from the above system can be 
written as 

d a d a (f> + R^E H + a 2 d a d a <J) + abd a d a X = 0, 

and 

b 2 d a d a X + abd a d a ct) = 0. 
These equations can be rewriten in the following form 



d a d a c 



0. 



and 



ad Q d a (j) + bd a d a X = 0. 
For static field configurations they become 



and 



dx 2 = ****** 
dx 2 dx 2 



0. 



(3) 
(4) 

(5) 
(6) 



According to the standard route to solitons pT| , the above equations of 
motion are the equations we have to solve to find soliton solutions. 

Before following this route, however, let us investigate the energy corre- 
sponding to static field configurations. Here we use the Lagrange density (Q) 
to write the energy as E = Eb + E', where Eb is 



E B = H[<f>(oo)\-H[<f>(-<x>)l 



and E' has the form 



E' = - 
2 



dx 



Tx~ H * 



/ d(p dx" 

a Tx +b Tx. 



2U- 



'dH* 



2-\ 



(7) 



(8) 



Here we recognize that the choice (0) for the potential lead us to the result 
that the energy corresponding to static field configurations can get to its 
lower bound, E B , and this happens when one sets 

dip dH 

dx dd> ' 



(9) 
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and 

a di + h di = {) - (10) 
As we can check straightforwardly, the above first-order equations (PD and 



(10) solve the equations of motion (y) and (y). These results are interesting 
since we can search for soliton solutions in coupled nonlinear (but now) first- 
order differential equations, and we already know that such solitons present 
minimum energy. 

Here we notice that xi x ) = X ~ ( a /^)0( a; ); with x constant, is the choice 
that solves the first-order equation (|T0|) . This is the solution we use for the 
X field. Furthermore, since our system is Lorentz- invariant, we can get time- 
dependent solutions by just boosting the static solutions. Hence, in order to 
mantain Lorentz invariance we must assume that x{ x ^) = X ~ (a/b)<p(x,t). 
In this case both 4>(x,t) and x{ x ^) satisfy the time- dependent equation of 
motion (|j). This is an interesting result, because now we can investigate 
classical or linear stability of the soliton solutions by just examining the 
equation for the (f) field. 

To make this investigation explicit, let us now write 

(f)(x, t) = (j)(x) + T](x, t), (11) 

where r](x,t) is the fluctuation about the classical solution <p(x). We substi- 
tute the above (f)(x, t) into the time-dependent equation of motion @ to get, 
up to first-order in 77, 

d 2 Tj d^Tj 

In this case we recall that the choice x( x , t) — X~( a /b)4>(x, t) is fully compat- 
ible with the investigation we are now dealing with. Furthermore, by assump- 
tion the classical solution <p(x) is static, and so the term H^ + H^H^ in the 
above equation is time- independent. Hence we can write rj(x,t) = rj(x)T(t), 
in order to separate space and time in the fluctuation. For this reason we 
can investigate classical or linear stability || by examining the following 
time-independent Schrodinger operator 

d 2 

82 = ~~^2 + H l<t> + H<l>H<t>4,<t>- (13) 



5 



On the other hand, we can use the first-order equation (|9|) to introduce 
the following operators 

Sf = ±^ + H H . (14) 

It is not hard to see that these first-order operators are adjoint of each other. 
Furthermore, they can be used to write the second-order Schrodinger opera- 
tor as S2 = Si Si . These facts ensure that S2 is positive semi-definite, and 
hence the static field <p(x) that solves the first-order equation (|9]) is classically 
or linearly stable. 

Before ending this section, let us now briefly comment on some topolog- 
ical issues. Since we are working with a bidimensional spacetime, we can 
introduce the following vector 

j a = e aP d p H(<f>), (15) 

where e af3 is the Levi-Civita tensor, which obeys e 00 = e 11 = and e 01 = 
— e 10 = 1. Hence we notice that d a j a = 0, and so there is a conserved quan- 
tity, the so-called topological charge, which can be written in the following 
form 

Q T = / dxf = #[0(oo)] - H[<j)(-oo)]. (16) 

J —00 

Here we notice that the current density ([15]) is defined |10| in terms of 



H = H((f)), in order to make the topological charge equal to the energy 
of the static field, as given in (|T6|). Evidently, other definitions for the cur- 
rent density can be used, but they make sense if and only if the potential 
has two or more absolute minima ]T1|. This is so because the field config- 



uration has to get to the absolute minima of the potential, asymptotically, 
to make the energy finite. Within this context, the definition we are using 
is very natural, because it uses H((fi), and H((fi) is just the function that 
defines the potential. For potentials that have only one absolute minimum, 
nontrivial field configurations must go to that value, asymptotically. Such 
solutions can not have topological charge, and this is the reason to name 
them nontopological solutions. 



3 Two-component hydrogen-bonded chains 

In this section we are concerned with mapping the class of systems introduced 
in the former section to hydrogen-bonded chains. Toward this goal, let us 
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first recognize that when we set a = b = 0, the above systems can be used to 
map the continuum version of one-component models for hydrogen-bonded 
chains. In this case the real scalar field <fi has to be identified with the proton 
degree of freedom in the unidimensional chain. 

To model two-component systems, we can naturaly use the scalar field 
<p to map the proton degree of freedom, in the proton sublattice. Within 
this context, we now take the scalar field \ to map the acoustic mode of the 
heavy-ion sublattice. Hence we notice that our class of systems can describe 
hydrogen-bonded systems with different on-site proton potentials, but with 
qualitatively the same linear coupling between the proton sublattice and the 
acoustic mode of the heavy-ion sublattice. As we have already shown, such 
a coupling leads to the solution x{x,t) = x — ip/a) (f>(x,t), where 4>(x,t) is 
obtained by boosting the static field 0(x) that solves the first-order equation 
(|9p, and hence the static second-order equation of motion. 

In order to illustrate the procedure, let us consider some explicit examples. 
Firstly we introduce 

F^ = Asin(0), (17) 
where A is a real parameter, with dimension of energy. Here we get 

H l = -Acos(0), (18) 

and so the potential has the form 

f/ 1 (0) = iA 2 sin 2 (0). (19) 

This is the sine-Gordon system, and the first-order equation (|9]) becomes 

^ = Asin(0), (20) 

and presents, for instance, the following soliton solutions 

4>(x) = ±2 arctanexp(Ax). (21) 

Here we notice that sign(A), which is unimportant to specify the system, 
gives the corresponding soliton and antisoliton solutions. See that the above 
solutions connect pairs of adjacent minima in the set of minima energy states 
= and <j)± = ±7r of the system. See also that there exist an infinity set 



7 



of solutions, which connect pairs of adjacent minima in the infinity set of 
minima energy states. As one knows, the sine-Gordon system is very rich, 
and there exist more solutions. 

As another system, let us choose 

H 24l = \{v 2 -<P 2 ), (22) 

where v is real and positive. In this case we have 

H 2 = A0 (V - l -<f) , (23) 
and the potential is given by 

U 2 {<P)= l -\\<f-v 2 f. (24) 
This is the 4 system, and the first-order equation ([)]) becomes 

! = P5) 

which presents the kink solution 

4>(x) = v tanh(At> x). (26) 

Here we also notice that sign(A), which is unimportant to specify the system, 
gives the kink (A > 0) and antikink (A < 0) solutions. See that these solutions 
connect the two minima energy states (0± = ±v ) of this system. 
As yet another system, let us choose 

tf30 = A0(u 2 -0 2 ). (27) 

In this case we have 

H«= ... , 

2 T V 2 



tfs^A^-V), (28) 



and the potential is given by 

1 



f/ 3 (0) = ^AW-^) 2 - (29) 
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This is the 6 system, and the first-order equation ([)]) becomes 



^ = A0(^-0 2 ), (30) 
It presents the kink solutions 

r l i 1/2 

0±(x) = ± j-t; 2 [l ±tanh(Aw 2 a;)]j . (31) 

Here we notice once again that sign(A), which is unimportant to specify the 
system, also gives all the solutions. See that these solutions connect pairs of 
adjacent minima in the set of three minima energy states (0 = 0,4>± = ±f ) 
of this system. 

The above two first examples are standard models for the on-site potential 
in hydrogen-bonded chains - see, for instance, Ref. ||. In particular, in || 
a more general potential was also considered. Here we refer to the following 
choice 

Tj x cos(0) - di . . 

4<p 1 - d\ 1 - d 2 [cos(0) - di] ' 1 ' 

where d\,d% are real and dimensionless parameters, with < d\ < 1. We 
notice that when d\ — di — 0, we get back to the sine-Gordon model already 
investigated. Yet, for <i 2 = we get to the double sine-Gordon model. 

Despite the complexity of H(4>), we can state on general grounds that 
for systems bellonging to the above class we have to deal with a first-order 
equation, instead of the usual second-order equation of motion. As we have 
already shown, we know that the soliton solutions are stable and have min- 
imum energy. The energy is obtained from (|^), and so we just need H(<p) 
and the asymptotic values of the field. After recalling that a topological 
soliton connects two adjacent minima of the potential, we immediately get 
the asymptotic values of the field, from which we write the energy easely. 
As an illustration, the above can be integrated to give if 4 , which is an 
important quantity. In the one-component model of Ref. || this plays the 
role of the general function Ak there introduced in Sect. II. E. Here, however, 
we have shown explicitly that the solitonic solutions have minimum energy, 
and are classically stable. 

The model fl27|) with the 6 potential is a new example. In this case 
there are the symmetric minimum (0 = 0), and two asymmetric minima 
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(4>± = =bu), and so this system seems to be richer then the 4 system, in 
which we have only two asymmetric minima, namely (0± = ±v). We have 
added this example to further illustrate the general procedure. 



4 Other two-component models 

The class of systems we have being investigating is specifyed by the potential 
@. However, we may find reasons to abandon this class of systems to con- 
sider other models. Such a situation has occurred, for instance, in Ref. [0], 
where an asymmetric </> 4 potential is used to find bell-shape soliton solutions. 
This is the case we want to consider in the following. Before doing this, how- 
ever, without loosing generality let us first shift the <fi field in the 4 potential 
given in the former section. Here we change <fi —>■ <fi ± v to get 

U{4>) = 2AV0 2 ± 2\ 2 v<f> 3 + ^A 2 4 . (33) 

Of course, this is still the potential for kink solitons we examined in the former 
section, shifted by the value we have just introduced. Here the solitonic 
solutions are given by 

<f>±(x) = T v + vtanh(Xvx). (34) 

To get to the case of an asymmetric 4 potential, let us now add to the 
above potential (|33|) the following contribution 



U'(<f>) = -2eAV0 2 , (35) 
Here e is real, positive, and dimensionless. In this case we obtain 

U{(j>) = 2(1 - e)AV0 2 ± 2AV 3 + ^A 2 4 , (36) 

and so it does not bellong to the class of systems we have introduced in 
Sect. 0. To search for solutions we have to deal with the equation of motion, 
which is, for static configuration, 

^ = 4(1 - e)AV0 ± Q\ 2 v<p 2 + 2A 2 3 . (37) 
ax 2 
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This equation can be integrated to give 



1 + V e cosn(2Ai;vl - es) 

Here we notice that the above bell-shape solitons are nontopological solu- 
tions, and so we do not have any good reason to believe on their stability - 
recall that the proof of classical stability already introduced does not work in 
this case. For this reason, stability of the above bell-shape soliton solutions 
has to be examined independently, but this is out of the scope of the present 
paper. We shall return to the issue concerning the presence of nontopological 
solitons in systems of coupled scalar fields in another paper. 



5 Ending comments 

In this paper we have shown explicitly that the class of systems introduced 
in Sect. || can be used to model hydrogen-bonded chains. As we have also 
shown in Sect. ^, this form of modelling two-component systems of hydrogen- 
bonded materials is specific, and works if and only if the two components 
couple linearly, in the form of a derivative coupling. 

The most interesting results we have found in the present route to soli- 
ton solutions in hydrogen-bonded chains can be summarized as follows. To 
search for solitonic solutions, it sufices to deal with first-order differential 
equations, which solve the corresponding equations of motion for static con- 
figurations. The solitonic solutions one may find present minimum energy 
and are classically or linearly stable. In addition, we have also introduced a 
way of accounting for topological properties of solitonic solutions. 

The above results encourage us to search for more general systems of cou- 
pled scalar fields, with the aim of perhaps finding classes of systems that lead 
to a unified way of investigating two-component models of hydrogen-bonded 
chains, without assuming linear coupling between the two components from 
the beginning. This and other related issues are presently under considera- 
tion. 

We would like to thank R. F. Ribeiro for some interesting discussions. DT 
is grateful to Conselho Nacional de Desenvolvimento Cientffico e Tecnologico, 
CNPq, Brazil, for a fellowship. 
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